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Abstract 

Silicon nanowires grown along the (lOO)-direction with a bulk Si core are studied with density 
functional calculations. Two surface reconstructions prevail after exploration of a large fraction 
of the phase space of nanowire reconstructions. Despite their energetical equivalence, one of the 
reconstructions is found to be strongly metallic while the other one is semi-metallic. This electronic- 
structure behavior is dictated by the particular surface states of each reconstruction. These results 
imply that doping is not required in order to obtain good conducting Si nanowires. 

PACS numbers: 73.22.-f, 81.07.Bc, 81.07.Lk 
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One-dimensional semiconductor nanostructures are attracting great interest for their po- 
tentially high impact in future molecular-electronics applications, such as nanoswitches and 

n 

nanocontacts [Ij. Silicon nanowires (SiNWs) appear to be an especially appealing choice, 
due to their ideal interface compatibility with conventional Si-based technology P,!^. More- 
over, measuring changes of their conductance allows to use SiNWs as real-time label-free 
sensors in different chemical environments. The detection of NH3 4] and of biological macro- 
5| has been reported, envisaging the possibility to reach the single-molecule de- 



molecules 
tection limit 



n 



Recently, SiNWs of diameters below 10 nm have been synthesized. Ma et al. 7] have 
reported extremely thin SiNWs grown along the (110) lattice direction, whose diameters 
range from 1.3 to 7 nm. Previously, Holmes et al. [8] obtained (100) and (110) SiNWs of 4 
to 5 nm and discussed the influence of lattice orientation. Wu et al. |9j have grown (110), 
(111) and (112) SiNWs down to 3 nm. Colemann et al. have reported SiNWs of 3 to 
5 nm diameter. 

Thus, a detailed understanding of thin SiNWs structure and of their mechanical and 
electrical properties as a result of the different growth conditions is required. While it has 
been extensively demonstrated that in H-terminated SiNWs quantum confinement induces 



a gap-broadening effect 
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13|, little is known about surface reconstruction of non- 



passivated wires and about their electronic structure. Theoretical studies of the structure 
of the thinnest possible SiNW Q, Q| have been previously published. However, recent 
experiments have provided convincing evidence that some SiNWs grow around a mono- 
crystalline bulk Si core l^.lla]. Systems fulfilling such a requirement have been studied by 
Ismail-Beigi and Arias |n| ((100) SiNWs) and by Zhao and Yakobson Q ((100) and (110) 
wires). Both studies show the paramount importance of facet edges in wires with diameters 
in the nanometer range. 

Thin nanowires pose fundamental problems with respect to their conduction properties. 
Recently, it has been discussed that doping may not be the advisable technique for tailoring 
nanodevice electrical properties, because of the expected statistical deviation of impurity 



concentration from one system to another 



taining the measured SiNWs conductance 



21 



3 



20|. Indeed, typical concentrations for at- 
may well mean that no donor / acceptor is 



actually present in nm-wide SiNWs. Some of these measurements show that conductance is 
larger than expected for doped SiNW 121| and the authors suggest that surface states may 
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be responsible for the measured conduction. 

In this Letter, we present a theoretical study of realistic (100) SiNWs, with a bulk Si 



core and a diameter of ~ 1.5 nm 
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18|, thoroughly exploring the phase space of SiNW 



reconstructions. We study different periodic cells, ranging from 57 to 912 atoms, so that pe- 
riodicity does not prevent us from finding the lowest-energy reconstruction. The mechanical 
properties of the SiNWs are studied by computing their Young modulus and Poisson ratio. 
Particularly, we focus our attention on the electronic properties of the wire as determined by 
the different possible lateral facet reconstructions. We give the first description of nanowire 
surface states as they evolve from the particular surface reconstruction of the studied SiNW. 
Depending on the reconstruction, the surface state can cross the Fermi level affecting the 
electrical characteristics of the nanowire. Hence, there is an intrinsic relation between the 
reconstruction of nanowire facets and its transport properties. This is a remarkable finding 
that shows that doping is not needed for obtaining conducting SiNWs. 

We have performed density-functional theory (DFT) calculations with both a numerical 
atomic orbital and a plane- wave 12J| basis set. We have used a double-C polarized basis 

n n 

set |23[ and a plane- wave energy cutoff of 20 Ry [2J|, with the Generalized Gradient Ap- 
proximation 25] for the exchange-correlation functional. We have studied wires in supercell 
geometry with a diameter of ~ 1.5 nm. The axis periodicity will restrain the number of 
possible reconstructions and for this reason we have considered different supercell sizes, an- 
alyzing SiNWs of 57, 114 and 171 atoms. The reciprocal space has been sampled with a 
converged grid of 12, 6 and 4 k-points respectively. The atomic positions have been relaxed 
until the maximum force was smaller than 0.04 eV/A. 

The faceting geometry adopted by the wire is given by thermodynamical considera- 
tions HQ , |l9|. On the one hand, the formation of {100} facets is favored over {110} 
facets, due the lower corresponding surface energy; on the other hand, facets with an even 
number of atoms can dimerize, lowering their energy, and are thus favored over facets with 
an odd number of atoms [2^. 

We have obtained two competing geometries for the {100} facets: a Ic reconstruction [see 
Fig. (^)] and a 2c reconstruction [see Fig. ^p)], being c the lattice vector along the axis of 
the unreconstructed wire. The Ic reconstruction has the same periodicity of bulk Si, presents 
a trough in the middle of the facet and turns out to be the most stable, with a cohesive energy 
of 3.990 eV/atom j^^j. This value is only 3 meV/atom lower than for the 2c reconstruction. 
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The small energy difference - probably within the accuracy of our calculations |28 . 
indicates that both reconstructions are likely to coexist. The trough in the Ic reconstruction 
is determined by the buckled dimer sequence, which presents a double chain of low atoms 
[labeled with a circle in Fig. ([Q)]; this rule is respected by only one of the two sides in the 2c 
reconstruction, while on the other side one every two dimers is flipped. In contrast to what 
has been proposed for the Si(lOO) surface js^, none of the two cases shows spin polarization. 

The two arrangements do not correspond to any infinite surface reconstruction, even 
though they loosely resemble the Si(lOO) c(4 x 2) and p(A x 2) (but not the p{2 x 2), for the 
Si(lOO) surface reconstructions see for example Ref. [221). The difference between surface and 
SiNW reconstructions stems from the lower coordination of the facet atoms. Between two 
adjacent dimers on the facet there is one single atom in the underneath layer, while between 
the corresponding dimers on the (100) surface there are two. This reduced coordination 
leads to a lateral shift of the dimers, increasing their packaging. 

In order to extensively explore other possible reconstructions, we have also performed 
non-orthogonal tight-binding (TB) j^^l^^ calculations, considering supercells with a lattice 
parameter along the wire axis up to 16c (912 atoms). The TB structural relaxations are in 
very good agreement with the DFT results and no new reconstruction was found. We have 
also taken into account other faceting arrangements prior to relaxations, but we confirm 
that the minimum energy configurations are obtained when {100} facets prevail 17, 18, loj ]. 
Second-neighbor empirical potentials have also been tested, but they have proven 

to be unsuitable to reproduce the complexity of these reconstructions. 

Given the small energy difference between the two reconstructions, we have checked if at 
finite temperatures one of the two phases prevailed more clearly. We have calculated the 



Helmholtz free energy F following a quasi-harmonic approach up to 300 K [37|. We have 
found that the difference between the two structure remains practically unchanged all over 
the temperature range analyzed. This is due to the fact that the vibrational contribution to 
F is determined by the integrated phonon density of states which hardly changes from one 
case to the other. 

If SiNWs are to be used as nanoswitches or for manipulation purposes, it is important 
that they have a certain stiffness that prevents them from collapsing as a result of mechanical 
tensions. We have studied the response to axial stress in (100) SiNWs within non-orthogonal 
TB, finding a Young's modulus of ~ 137 GPa (~ 195 GPa for bulk Si) and a Poisson ratio of 
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~ 0.35 j38|. These values confirm the intuition that their bulk core make SiNWs mechanically 
stable. 

We have calculated the band structure diagrams corresponding to the two different recon- 
structions. Surprisingly, the electronic structures turn out to be utterly different, especially 
if compared to their rather similar geometries. The Ic-reconstruction has four clearly metal- 
lic states [see Fig. ^^)], while the 2c wire shows a semi-metallic behavior with a single 
Fermi-level crossing at the zone boundary [see Fig. (jib)]- The metallic character of the Ic- 
wire is robust against failures of DFT. A scissor operator correction would rigidly shift some 
of the electronic bands. In order to open a gap, both of the degenerate bands (ii) and (iii) 
in Fig. (j2K) should become fully occupied or fully empty because both bands have the same 
symmetry (see Fig. (jH))). Given the band structure about the Fermi energy, it is not possible 
to find shifts of other bands to prevent this shift from violating electron counting. Thus the 
Ic-wire has metallic character. However, the limitations of DFT approximations to give the 
correct band gap of bulk silicon only permit to claim that the 2c-wire is a semi-metal of a 
small gap material. 

The two reconstructions evince two competing mechanisms ruling the surface energetics 
in connection with the electronic structure of the SiNWs. There is a trade off between 
recovering the bulk-like tetrahedral coordination on one side and the full delocalization of a 
genuine Bloch state along the wire axis on the other side. The first of these two mechanisms 
favors the flipping of one of the dimers which is characteristic of the 2c-reconstruction. The 
larger coordination of the protruding Si atom in the flipped dimer leads to angles of ~ 100°, 
while the corresponding angles in the Ic-reconstruction are of ~ 90°. This energy gain is 
compensated by the shift of the electronic structure as shown in Fig. (0). This shift is due 
to the reduction of symmetry by the flipped dimer. The flipped dimer interrupts the surface 
Bloch states, hence the overlap between dangling bonds is reduced in the 2c-reconstruction 
and the energy of the surface states increases. Electronic localization induces an increase 
of energy by breaking the one-dimensional Bloch state, but on the other hand makes the 
system lower its energy locally by increasing the bonding of one of the dimers. 

The bands of the Ic-reconstruction can be identified with the 2c-reconstruction ones. 
Fig. (0). The (ii) band of Fig. (j2^) shows similar symmetry as the (2) of Fig. (03) (see 
Fig. 121 below). An adiabatic calculation, slowly flipping the dimer yields this one-to-one 
correspondence. In the same manner, bands {iii) and {iv) of Fig. (j2K) can be identified with 
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(4) and (3) of Fig. (jJb) respectively. The effect on tlie electronic structure of flipping one 
dimer is to localize, leading to higher-energy and flatter bands, Fig. ^p). 

The electronic structure about the Fermi energy is given by the surface states originating 
in the dimer dangling bonds. This is shown in Fig. (jH)). Iso-surfaces of wave- function 
amplitude have been plotted on the SiNW atomic structure, showing that Bloch states 
originate on the facet dangling bonds. Figure Q also shows the correspondence of the 
wave functions (a)-ii and (b)-2, as explained above. Surface states are thus present in both 
reconstructions. 

Conductance through thin SiNWs remains a difficult issue The way SiNWs bind to 
metal electrodes is complex and crucial for determining the actual conductance of SiNW- 
based devices. In the absence of defects and in the case of excellent electrode contacts, 
the strongly metallic SiNW (Ic-reconstructed) will present a maximal conductance of four 
quanta, equivalent to a resistance of only ~ 3kQ per SiNW. In the same conditions the 
semi-metallic SiNW (2c-reconstructed) exhibits a minimum resistance of ~ 13kQ per SiNW. 
These results show that SiNWs can be good conductors, and their actual electrical properties 
will be strongly dependent on the growth conditions. Experimental studies 0| show that 
despite diffusion from the metallic contacts and/or surface contamination the measured 
conductance of thin SiNWs after annealing is much higher than the one expected from 
doping the SiNWs. This finding is in qualitative agreement with the above surface-state 
driven conduction. 

Any perturbation on the surface states will drastically affect the conduction properties of 
the wires. These wires will be extremely sensitive to the chemical environment and therefore 

a.e,oodca.d,date.fo..,no>ecu>a..detect,o4Ba. 

In conclusion, we have performed DFT calculations allowing for different reconstructions 

of (100) SiNWs with a bulk core. The SiNWs are classified according to their {100}-facet re- 
constructions. No direct correspondence is found with infinite surface reconstructions, hence 
the two minimum-energy possibilities are classified according to the size of their periodicity 
along the wire axis: Ic for the one with the shortest period and 2c for the longest one. 
The bulk-like core of the wire confers them with a Young modulus close to the bulk Si one; 
on the other hand, the Poisson ratio indicates good lateral elastic properties of these wires. 
Despite the closeness in energy of the two reconstructions, the Ic one shows a strong metallic 
character while the 2c one is semi-metallic. We have rationalized these findings in terms 
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of the gain in energy due to the formation of surface state versus tetrahedral-hke angles in 
the surface dimers. Surface states are probably ubiquitous in a large family of SiNWs. The 
wires studied here present good conduction properties thanks to the reconstruction-induced 
surface states. Hence, the use of thin SiNWs as conductors without doping is possible. 
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FIG. 1: {100}-facet reconstructions, (a) Ic and (b) 2c dimerization of the surface dangling bonds. 
Lattice vectors of the facet unit cell are sketched to underline the different periodicity along (100) 
axis (c is the bulk lattice parameter). 
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7c/2c 
(a) 
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(b) 



FIG. 2: Band structure diagrams, (a) The Ic reconstruction has four metalHc states [(i) to (iv)] 
clearly crossing the Fermi level; (b) the 2c reconstruction has rather a single semi-metallic state. 
For comparison, both diagrams refer to the 2c unit cell. 
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(a)-iv 



(b)-2 



FIG. 3: Cross-section view of the wave functions at T point of (a) the four metahic states of the Ic 
reconstruction [(i) to (iv)] and (b) the semi-metalhc state of the 2c reconstruction. Labehng follows 
the scheme of Fig. These wave functions give a false impression of C2v symmetry. The pseudo 
C2v symmetry explains why bands {ii) and (iii) are almost degenerate, (surface states (a) — ii and 
(a) — iii in the figure) . 
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